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Abstract 



oo 

o 

, The perturbative effective potential V in the massless A<^'* model with a global 0{N) symmetry is uniquely 

determined to all orders by the renormalization group functions alone when the Coleman- Weinberg renormalization 

d^v I 



^ , condition | = A is used, where fi represents the renormalization scale. Systematic methods are developed 

pL^ I to express the n-loop effective potential in the Coleman- Weinberg scheme in terms of the known n-loop minimal 

subtraction (MS) renormalization-group functions. Moreover, it also proves possible to sum the leading- and 
' subsequent-to-leading-logarithm contributions to V. An essential element of this analysis is a conversion of the 

renormalization group functions in the Coleman- Weinberg scheme to the renormalization group functions in the 
MS scheme. As an example, the explicit five-loop effective potential is obtained from the known five-loop MS 
' renormalization group functions and we explicitly sum the leading logarithm (LL), next-to- leading (NLL), and 

further subleading- logarithm contributions to V. Extensions of these results to massless scalar QED are also 
^ • presented. Because massless scalar QED has two couplings, conversion of the RG functions from the MS scheme 

' to the CW scheme requires the use of multi-scale renormalization group methods. 

^ ■ 1 Introduction 

' The effective potential V in the massless A^"* model with a global 0{N) symmetry has received considerable attention 
\^ , [U [21 [3l m [5] because of its connection to the scalar-field theory projection of the Standard Model for iV = 4. A 
C ' variety of renormalization schemes have been employed in computing the effective potential, among them minimal 
subtraction (MS) [S] and imposing the Coleman- Weinberg (CW) condition [TJ [31 [7] 



= 24A , (1) 

0=M 



o ■ 

where ^ is the renormalization scale. 

In this paper we develop iterative techniques that uniquely determine leading-logarithm and subsequent-to- 
leading-logarithm expansions of the effective potential in the CW scheme for 0(iV)-symmetric massless Xcf)'^ theory 
' and for massless scalar QED with a interaction. As discussed below, the renormalization-group (RG) functions 
in the CW scheme differ from the known minimal-subtraction (MS) scheme RG functions |7j, resulting in non-trivial 
effects of this scheme conversion beginning at two-loop order. Although it has been known for some time that the 
effective potential is in principle determined by the RG equation [l] , two-loop calculations have either failed to make 
the necessary scheme conversion ^ i9j or have been done explicitly without using RG methods ^ . 

In Section [5] we explicitly construct the effective potential V for the X4>'^ model, not only up to five loop order, 
but also give all the N^LL (next-to, next-to, next-to, next-to leading-logarithm) contributions to V without any 
explicit evaluation of any diagrams, simply by applying the RG equation in conjunction with the CW renormalization 
scheme, thereby realizing the result of Ref. [T]13 

In Section [3] we extend our analysis to massless scalar QED with a cj)'^ interaction, a theory which contains two 
couplings. The results are quite similar to the A(/)^ scenario; iterative methods are developed to determine the scalar- 
field effective potential in terms of the RG functions in the CW scheme. The presence of multiple couplings requires 
the use of multi-scale RG methods [101 [H] to convert the RG coefficients to the CW scheme. 
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^The relation between V and the RG equation appearing in Ref. 1 is further analyzed in Ref. |12| where the importance of fixing the 
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2 Massless 0(A^)-Symmetric Theory 

In 0(A^)-synimetric massless theory, the effective potential in the CW scheme takes the form 



oo n 



V{X, 0, ^) = ^ ^ A"+ir„™L™</.4 (2) 

n— m— 

when computed in perturbation theory, where we have defined 

L^log(^). (3) 

Since the renormalization scale ^ is unphysical, changes in /i must be compensated for by changes in A and (f>; this 
leads to the renormalization group equation 

o o Q \ 

M^+/?(A)7TT+7(A)0Tr7U^(A,0,M)=O (4) 



9/i dX 



where 



/3(A)=/i— = ^6feA'= (5) 

fe=2 

The renormalization group equation (j4]) can be used [3 I13j to sequentially sum the logarithms appearing in 
Eq. 12]); that is if we rewrite Eq. ([2]) as 

oo 

where 

oc 

Sn{XL) — T„+„i„, (AL)™ , (8) 

then the functions SniO determined by Eq. ^ provided we impose the boundary condition 5'„(0) = T„o- More 
explicitly, we write Eqs. (UHZ]) as 

oo 

= ^ [-1 + (.giA + g^X^ + ...)] iX^+'S^iO + {b^X^ + b,X^ + . . .) X'^ [{k + 1)3^(0 + ^S^O] 

k=o (9) 
+ 4(giA + 52A2 + ...) A'=+i5fe(0 , 



so that to order A'^, and in general to order A" we respectively find that 

So = (10) 



(-2 + 4 + (^2 + 4.gi) 



= 



(-2 + &2O ^ + {n + l)b2+ 451 



In general, if Eq. Q is satisfied at order A"+^ , 5„ (^) satisfies the differential equation (fTT|) whose solution requires 
that we know S'm(C) ('t* = 0, 1, . . . rt— 1), fe„i (m = 2, . . . n + 2), t;™ (m =1, . . . n+ 1), and the boundary condition 
^^(O) — Tno- In other words, Sn{0 is governed by coupled differential equations which depend upon the n + 1 loop 
RG coefficients. It is important to note that the RG equation by itself does not determine the boundary conditions 
T„o; these will be seen to be determined by the CW renormalization condition. 
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The solutions for 5*0 and Si are 

SoiO = — (12) 
w 



4.927oo , 4gr27oo + b2Tio b^Ti 



S"! = 7 + T 2 T (13) 

O2W 02W^ 02W^ 



where w = 1 ~ and gi = for the A(/)^ model As suggested by p2)) and (fT3)) . the expUcit solutions to Eqs. (fTOl 
and (jlip take the form of polynomials in w and logw: 



i=l j=0 



Expressions for the coefficients a"j up to n = 2 are given in Appendix A in terms of a recursion relation. Appendix 
A also demonstrates that a partial summation of the recursion relation is possible. 

If V is computed so that Eq. ([T]) is satisfied, then upon substituting Eq. ([7]) into Eq. ^ , we find that 



24A = 



(15) 



(We note that L = if ^ = /i.) It follows from Eq. ([T5|) that for the CW renormalization scheme 



^o(O) = Too - 1 (16) 

25 

= 1005^(0) + 245i(0) = 5O62T00 + 24Tio =^ = -— ^2 (17) 

= 140S'^'(0) + 1005^(0) + 245-2(0) = 12r2o + 5O62T10 + (356^ + 2563 + IOO52) Too (18) 
=80S'^"(0) + 140S';'(0) + 1005^(0) + 24^3(0) 

(6O53 + 17562^3 + 5O64 + 6105232 + 2OO53) Too (^g^ 

2106^ + IOO63 - 48^ + 200^2 - 48f^ + ^ [64 + 53] ) Tio + I5O62T20 + 24T30 



62 62 62 

etc.] in general by having Eq. pS]) satisfied at each order in A we end up with 



16^5fc(0) + 80— 5fc+i(0) + 140-— 5fc+2(0) + 100— 5fe+3(0) + 245fe+4(0) = (A: = 0,l,2...) . (20) 



Consequently, the (rt + 1 loop) boundary condition ^^(O) = T„o is determined iteratively by lower-order results via 
Eq. ((20)) : that is, once Sk{£,) ... Sk+ziC) are known, S'fe4.4(0) = Tfe+40 is fixed by Eq. ((20l) . Hence V is determined 
entirely by the renormalization group functions /3(A), 7(A) when employing the CW renormalization condition of 
Eq. ([T]). It is not apparent how T„o can be determined in any other scheme except by relating that scheme to the 
CW scheme. 

However, the renormalization group functions are generally computed in the MS renormalization scheme; being 
given up to five loop order in Ref. |14j . To relate the renormalization group functions in these two schemes, we note, 
following Ref. [7 , that in the MS scheme, the form of the effective potential is much the same as that of Eq. ([2]), 

00 n 

V{X, ^,fi) = Y.Y. A"+'T„„L"V" (21) 

n— m— 

where now 

~L = log (^) . (22) 

Upon comparing Eqs. ([3]) and ((22)l . the mass scale /i^ in the MS scheme can be related to the mass scale /x^ in the 
scheme of Eq. ([1]) by 

il = X^fj, (23) 
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Consequently, ^ = - ^^/3(A) where /3(A) = jl^, 07(A) = /i^ and thus [7] 

/3(A) = (24) 

2A 

relates the renormalization group function in the two schemes. Knowing /3(A) and 7(A) in the MS renormalization 
scheme thus determines /3(A) and 7(A) in the CW renormalization scheme [Eq. ([T])] and hence V can be determined 
entirely from /3(A) and 7(A). In particular, if /3(A) = b2X^ + h^X^ + . . ., 7(A) = gi\ + g2\^ + . . ., then Eqs. (pi]) 
and ()25p can be expanded to convert the five- loop MS-scheme renormalization group functions of [14j . The explicit 
results to three-loop order (the first order at which the conversion of the anomalous dimension is non-trivial) are 

62 = &2 (26) 

63 = &3 + \bl (27) 

64 = 64 + + (28) 

51 - 31 = (29) 
92 = 32 (30) 

53 = 33 + ^^252 (31) 

The results up to five loop order (or indeed any desired order) are easily obtained. 

The MS-scheme RG coefficients for the 0{N) version of the At/)'' model are known to five-loop order [14]; to 
establish our conventions their values to three-loop order are 

i3^-«^ (33, 

- 33iV2 ^ 922/V + 2960 96(5/V -f 22)C(3) , , 

^4 = (34) 

31 = (35) 

/V + 2 

h = (36) 

_ iN + 2)iN + 8) 

128^^ ■ ^^^> 

The sum of all leading-logarithm (LL) and next-to-leading-logarithm (NLL) contributions to V is given by 

Vll+nll = A [So{XL) + XSiiXL)] <j>^ (38) 

where 5*0 and Si are completely determined by Eqs. (HH), US]), 1^, Wc 
can then recover the complete two-loop CW-scheme result for V computed explicitly in Ref. [3] upon expanding 
Sq and 5*1 to terms quadratic in L = log ^ and calculating T2Q via (jlSp . It should be noted that the two- loop 
CW-scheme result ^ does not satisfy the RG equation with MS coefficients; if one applies solves the RG equation 
with MS-scheme RG functions it would disagree with the explicit two-loop calculation. 

As noted earlier, S4 requires knowledge of the renormalization-group functions to five-loop order, and hence ^4 
is the highest-order term in the expansion ([7]) that can be determined with current knowledge of the renormalization 
group functions in the massless 0{N) theory [TJ. The solutions 5*0, Si, S2, S^ and S4 contain the boundary-condition 
coefficients T„o for n < 4. These five coefficients are determined by the set of five equations pBjl ^ P^ and with 
A; = 0. Once these coefficients are determined, the N*LL expression for the effective potential 

4 

^w4ii = ^A"+i5„(AL)04 (39) 

n=0 
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is expanded up to fifth order in L to obtain the five- loop perturbative coefficients Tnm for < n < 5 and < to < 5. 
The remaining five- loop coefficient T^q = £'5(0) is determined via ((20)) with k = 1. We present our determination 
of the explicit values for the perturbative coefficients to three-loop order (one order higher than the CW scheme 
calculation of Rcf . [3] ) 



Too 
Tio 

720 

Tbo 

T31 

732 



1 



25(iV + 8) 



, Til 



247r2 ' 47r2 
5 (177V2 -I- 347iV -I- 1418) 



T21 



lliV2 + 2067V-l-836 _{N + %y 
, J 22 



727r4 ' 247r4 ' IGtt* 

5 [7847V3 _^ 26305iV2 + [251338 -I- 7200C(3)] N + 694032 + 31680C(3)] 

23047r6 

296^3 _^ 9425Ar2 ^ [37242 -I- 1440C(3)] N + 239376 -f 6336C(3) 



3847r6 



[N + 8) (l0iV2 + 2097V 
64^1^ 



, ^33 — 



{N + Sf 
647r6 



(40) 
(41) 

(42) 

(43) 

(44) 

(45) 



The results (HUt - P^ are in agreement with the explicit two-loop calculation [5]. As mentioned earlier, the corre- 
sponding expression 52 (^) used to obtain these coefficients is given in Appendix A. Although the analytic expressions 
for the remaining coefficients to five-loop order are too lengthy to be presented, we give their numerical values for 
= 1 (simple scalar field theory) and = 4 (the scalar-field theory projection of the Standard Model) in Tables [1] 
and El 





T40 


Til 


Ti2 


^43 


T44 


N = 1 


5.218 


-2.277 


0.4457 


-4.865 X 10^2 


2.701 X 10--' 


TV = 4 


14.59 


-6.477 


1.304 


-0.1475 


8.537 X 10-'^ 



Table 1: Four-loop perturbative coefficients for = 1 and A^ = 4. 





750 


^51 


T52 


7^53 


754 


Trr 

^ 55 


N = 1 


-14.06 


6.407 


-1.362 


0.1744 


-1.404 X 10^2 


6.158 X 10-* 


A^ = 4 


-50.09 


23.20 


-5.065 


0.6709 


-5.631 X 10-^ 


2.595 X 10-^ 



Table 2: Five-loop perturbative coefficients for A^ = 1 and A^ = 4. 

In the next Section, we examine how the methods developed for massless scalar field theory can be extended to 
massless scalar electrodynamics, a theory with multiple couplings. 



3 Massless Scalar Electrodynamics 

Massless scalar quantum electrodynamics (MSQED) has the Lagrangian 



(46) 



The effective potential V{(l)) in this model can be computed perturbatively in a variety of ways [TJ [H [31 IH [H] . The 
effective action for MSQED has the expansion 



In the CW renormalization condition [I] 

d4F(0) 



= 5, 



(47) 

(48) 
(49) 
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one finds a perturbative expansion for ¥{(/)) of tlie form 

(oo n-\-k oo \ 
E E E Tn+k^r^r^kg^+'-'-a^L" 0^ ^ (5O) 
n=l r=0 k=0 ) 

L^log(^), (51) 

where a = e^, and ^ is the RG scale appearing in Eq. (|48p. The effective potential satisfies the RG equation 

= = (,| + /^^l + r 1^ + a, 0, ,)^BV. (52) 

Here, /?^,/?" and 7 are the renormalization group functions 

,00 n 
,00 n 

/3"(5,«)=A^^ = E^» ' /3^=^6^,,,5"-V (54) 

7(5, ") = = E ^" ' '^''^ = E 7n-r,rff""''a'' (55) 

In the previous Section it was shown that the effective potential in an 0(Af)-synimetric massless At/i^ theory is 
uniquely determined by the RG equation in the CW scheme. We now extend this analysis to deal with the situation 
occurring in MSQED where two couplings g and a occur. As shown below, there are two crucial distinctions between 
the single- and multiple-coupling situations. First, the coupled ordinary differential equations (fTTj) get replaced by 
coupled partial differential equations. Secondly, the conversion of the RG functions from the MS scheme to the CW 
scheme requires use of multi-scale RG methods (TUl [TT] . 

We now proceed to shown how together Eqs. ([^5]) and ([5^ again determine V without the calculation of additional 
Feynman diagrams. We first define 

n 

p^(g, a)^Y. rn-r.r,fc5""''«'' > fc + 1) (56) 
As a result, Eq. (|50p can be written 

oo n — 1 

y(g,a,0,/i) = ^^pJ;(5,a)iV\ (57) 

n=l fc=0 

with the contributions 

oo 

^LL-E^'fe+i^''^"' (58) 

oo 

Vnll = Y.v\^'^^H^ ■■■ (59) 

oo 

^iv^iL = E^'tp+ii'</'' (60) 



fc=0 



giving the leading-log, next-to-leading-log etc. corrections to V . 

First, we find that substitution of (l57l) into (l48l) results in the condition 



24p° + lOOpi 280p^ + \mvl -t- 384p4 = g (61) 
upon noting that if </> = /i, then L = 0. Recalling that = if n < fc + 1, we see that if n = 1, Eq. ((6T|) leads to 

= ^ , (62) 
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the tree level result. 

We now substitute Eqs. ([531) -([551), dST]) into Eq. jSH) to obtain 



EE 

n=l fc=0 



-2fcp«L 



k T fe-1 



E 

m=2 



' dg 



da 



m — 1 



</.4 = 



(63) 



Since p^, P^, (3" and 7„ are all polynomials of degree n in g and a, we can obtain coupled first-order partial differential 
equations that express each of the in terms of /3" and 7„. This is done by requiring that Eq. ([63]) be satisfied 
order- by-order in L*"' and in X]r=o '^'"S""'^'^'^' These equations are first order partial differential equations whose 
boundary conditions are provided by Eq. (|6f p . 

For example, at second order in the couplings and order zero in L, Eq. ([SS]) leads to 



M + l3l^+P2^+ilip"i-0- 



dg 



dc 



(64) 



This fixes p\, since is given by Eq. 
that 



Now at third order in the couplings and order L, we see from Eq. 



.2 , fl9M+/3."M 

dg ^ da 



+ 471P2 = • 



(65) 



Since P2 has been determined, Eq. (|65p serves to fix p^. In general, at order (n + 2) in the couplings, and order n in 
L, we find that 



Pn+2 — 



1 



2(n-f 1) 



d 



d 



/32V+/32T^+47: 



dg 



da 



Pn+l ■ 



(66) 



so that all of the contributions to Vll in Eq. (|58p can be calculated provided we use the expression for pi given by 
Eq. (Ha). 

Since we now know P2 , we can set n = 2 in Eq. (|6ip , leading to 



P2 



25 1 



6 



■P2- 



(67) 



This serves to start the sequence P2,pl ■ ■ ■ etc. Upon looking at terms in Eq. that are of order n -I- 3 in the 
couplings and n in L we have the recursion relation 



Pn+2 - 7lP„+lj - TT- 



1 

2n 



d 



d 



^'^d-g^^^d-a 



Pn+l 



/33^+472)pr 



(68) 



■P 



P+i 



which finally fixes all contributions to Vnll in Eq. in terms of /3|, /?2 , (3^, 71 and 72. 

It is evident that Vatpll involves the polynomials Pn+p+i- From Eq. (ICTj) . Pp+i can be found once Pp+i 
have been computed in the course of determining Vjvp-ill . . . Vll- Having fixed Ppj^i in this way, all subsequent con- 
tributions to Vj\[PLL are determined in terms of the polynomials Pn+p ■ ■ ■ Pn+i ''^sll as /3f . . . /Sp+i , /Sf ■ ■ • /^p+i , 7i ■ ■ • 7pi 
by considering those terms in Eq. (I63[) that are of order n + p + 2 in the coupling, and order n in L. The effective 
potential for massless scalar QED is therefore completely determined by the RG functions in the CW renormalization 
scheme. In Appendix B, the sums appearing in Eqs. ([55]) . (|55|) for Vll and Vnll are evaluated in closed form using 
a variant of the method of characteristics. 

As in the case of 0(Af)-symmetric massless Ac/)'* theory, it is necessary to convert the RG functions from the MS 
scheme to the CW scheme. However, because there are two logarithms log {gcfp' /p?^ and log [acj)'^ / fi^) appearing in 
the MS-scheme perturbative expansion, a simple rescaling of the renormalization scale fl as in (|23p cannot convert 
these two logarithms into the single logarithm (|5ip. 

The presence of multiple incompatible logarithmic scales is known to cause difficulties when attempting to solve 
the RG equation in other applications. To circumvent these problems, the concept of multiple renormalization scales, 
one scale for each appearance of the traditional MS RG scale in the Lagrangian, was first considered in HD] . This 
method was refined in by associating a renormalization scale with each kinetic term in the Lagrangian, which 
in the case of massless scalar QED, will introduce two renormalization scales {Kg and Ka) resulting in a MS-scheme 
perturbation series for the effective potential containing two logarithms: 



= log ( 



(69) 
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With multiple renormalization scales, there will exist MS-scheme RG equations and RG functions associated with 
each scale [11] 



dV 



d 



d_ 



where 



(70) 
(71) 

(72) 
(73) 
(74) 

As outlined in Ref. [Tl], these multi-scale MS-scheme RG functions can be obtained from the 1/e poles in the 
(multi-scale) renormalization constants. These multi-scale RG functions can also be determined by reconstructing 
the effective potential in the MS scheme from the MS RG functions and the logarithm-free parts of V; at this stage 
the renormalization scale jl can be split into Kg and allowing for a determination of Eqs. (|72p - (|74p through the 
requirement that V be independent of both Hg and Ha along the lines of Ref. [16| . Furthermore, in the limit when the 
two scales coincide {Kg — Ka ~ jl), the multi-scale MS RG functions are related to the single-scale MS RG functions 
/3s, /3" and 7 via [H] 





dg 


/3^ = 


dg 

dKa 


4" 


da 

^'''d^g' 


/3a = 


da 

dKa 




^-K ^ 

'dKg 




d(t> 

— I^a 

dKa 



(75) 



The RG functions ((72|) -([74 | must also be consistent with the integrability condition [Di,D2] V = 0; this constraint 
combined with the boundary condition (|75p may also be used to determine the multi-scale RG functions |11| . 
It is now evident that the rescalings 



(76) 



will convert the MS-scheme multi-scale logarithms ([M]) into the CW-scheme logarithm (I5T|) . thereby enabling scheme 
conversion. The RG functions in the CW scheme can then be obtained from (17611 combined with 



_ _ dKg dg dKa dg 

P ~ i^~r' ~ — "a 1" 1^' 



dg_ 
dfi 



dn dKg dfi dKo 



da dKg da dKa da 



^ ^ ^ dfj, dKg ^ ^ dji dK 



dfi 

d(t> 
dfjL 



d(j) dK„ d4> dKa d(j) 
dfjL dKg dfj, dKo 



to obtain 



/3" = /3" 



1 + ^ 
25 

2g_ 



■P'a 
-Pa 



l = lg 



1 + 



2g 



+ 7o 



14- 



/3" 
2^ 



(77) 
(78) 
(79) 

(80) 
(81) 
(82) 



The above equations can be solved perturbatively for the coefficients of the CW-scheme RG functions in terms of the 
multi-scale MS-scheme RG functions!! As expected, to lowest order one finds that the CW-scheme and MS-scheme 
RG coefficients coincide so that the effects of scheme conversion enter at two-loop level. 



•^It can be verified tliat the scalar field theory scheme conversion results I I24I I— I I25I I are obtained from the a — > limit of the MSQED 
results Il80|l - ll82| l . In this limit, /3| — > /3, 7g — > 7, and all other multi-scale MS RG functions become zero. Inversion of the resulting 

expressions /39 = /3 [1 + /39 /{2g)\ and 7 = 7 [1 -(- /JS /(2g)\ lead to pa =j3/\l- /3/{2g)] and 7 = 7/ [l - /3/(2g)] consistent with Eqs. ll24t - 

| [25t . 
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4 Conclusions 

In summary, we have developed iterative techniques that uniquely determine, in terms of MS RG functions, leading- 
logarithm and subsequent-to-leading-logarithm expansions of the effective potential in the CW scheme for massless 
Ac/)** scalar field theory with a global 0{N) symmetry. In these techniques, the N^LL expression is governed by a 
coupled set of first-order ordinary differential equations containing the p+1 loop RG coefficients, and the boundary 
conditions for this system are determined by the CW renormalization condition. In this approach, it is essential 
to convert the RG functions from the MS scheme (in which they are known to five-loop order in 0(A'^)-symmetric 
massless scalar field theory) to the CW-scheme. 

The methods developed for the scalar field theory one-coupling case have been extended to massless scalar 
QED. The presence of two couplings does not change the essential features of the analysis; instead of coupled 
ordinary differential equations the N^LL expansions are determined by systems of first-order partial differential 
equations resulting from the RG equation and algebraic equations arising from the CW renormalization condition. 
Similarly, conversion of the RG functions from the MS scheme to the CW scheme in massless scalar QED is also more 
elaborate, and requires the use of multi-scale renormalization group methods. Although multi-scale RG techniques 
are not widely known, the necessary multi-scale RG functions can either be calculated directly by introducing a 
renormalization scale for each kinetic term (and hence propagator) in the theory and exploiting the usual relation 
between the RG functions and the 1/e terms in the renormalization constants, or they may be reconstructed from 
the single-scale MS-scheme RG functions in conjunction with integrability conditions related to the commutator of 
the RG operator associated with each renormalization scale [101 E] ■ 

We would like to extend our methods to computing the effective potential when the mass of the field (f> is non- 
zero. In particular, our analysis may allow us to correct the two-loop renormalization-group analysis of the Standard 
Model effective potential Ref. [5] which is in disagreement with the explicit two- loop calculation J/Tj. It would also 
be interesting to see if the effective potential in the MS renormalization scheme could be determined uniquely by the 
renormalization group functions. 
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The differential equation pT|) establishes a recursive relation for the Sn{0 resulting in solutions of the form 



Appendix A: Massless Theory 




(83) 



where w = 1 — ^ — XL, and L = log(w). 

In terms of this notation, the solution for 5*0 is 



So{£.) = I o"? = ^2 Too = h , 
02W 



(84) 



and the solution for Si is 




(85) 



where 



CTi^o = -4g27oo , crio = ^^Tw + 452^00 , (^1,1 = -hzToo . (86) 

For the higher-order S'„, recursive expressions for afj provide the most compact form. For example, the solution 
for 5*2 is: 




(87) 
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where 

2 



' - -^(^3 + 4g2)ai,o - 2roo53 , (88) 



^2 
^2.0 



- {halo + 452Cr2,0 + (^3 - 452)ct2,i + (^4 + 62<?2)Too} , (89) 

-^92^1,1 , (90) 

^3% = ^ {(3^3 + 4g2)^T!,o + 832CT2\o + (2^3 - ^92)(tIi + (264 + + 26232)^00 + 262r2o} , (91) 

^li = -^3(2^2% - al^) , (92) 

-63^1 . (93) 



j2 

-'3,2 



The differential equation pT|) combined with the form of the solution ([55)1 can be used to obtain a set of recursion 
relations for the coefficients cr"^ . One finds that this procedure yields 

n-l 

=b2{j + 1)<, + 1 +{{n~i + 1)62 + 45i}ct,I!j + ^ [6„+2-m(j + IX^ + I + (« - l)(62.g„-m + 6„+2-m)CTrij , . 
- + l){b2g„-rn + 6n+2-m)0-" ij + i + {■ign+l-m + {m - 1 + l)&„+2-m)o-I"j] 

We note that af^- = if i > ?i + 1, j > i — 1, z < or j < 0. The coefficients for ^3 and ^4 can be extracted from 
the recursion relation ([M)) as needed to determine the T„m given in Tables [T] and O It is immediately apparent that 
if i = n + 1 and j = n in Eq. (|94p , then 

45i<Vi.„ = (95) 

and so for consistency gi = 0, as is already known from explicit calculation. If now in Eq. (|94p . we set i = n + I, it 
follows that 



where p = — ^. Considering values of i less than n + 1 results in a recursion relation that requires knowing 64, 172, 
etc. 

For j = n — 1 , it follows from that 

f^n+l.ri = PCr;';;^__;^ (97) 

so that 

= (98) 

where crj q — ^2- As a result, in the expansion of there is a contribution 

1 ^ T n 

^^=5:E^"^^-"Vn^ (99) 

n— 

which is a geometric series whose sum is 

Vi = 4 uT (100) 

4!^„ + ^logH 

Consequently, the sum of the contributions that are of the highest order in L and ^ at the N"LL order of 
perturbation theory gives rise to a singularity in V appearing, not when w — 0, but rather when w + log(w) = 0. 
If now j = n — 2 in , we find that 



n 



Cr;j+l_„_l = p ( ^^-7Y'^"."-2 ~ '^n,n-l ] (101) 



which implies 



= -np" + l + + + np-'al, (102) 
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where trl g is given above. Upon expressing 



11 1 V- 

- + - + ... + - = hm > 

2 3 n x^i^ 



k=2 



^dy ^ lim 



y - y 



■dy 



we now find that V now has a contribution 
which then becomes 



^ rn—l 

= y A"+ii— hm 

n=0 



-np 



y-y 



■dy + np^ V2_o 



(103) 



(104) 



(105) 



/o 1 - y 

Since X^fe^o kx^^^ = X]fc°=o ~ (i-x)'^ ' ^^"^ ^'^'^ integral in the above can be evaluated in turn, leading to 

log(l - B) 



Vii 



w [w + ^ log w 



pi 



B 



(106) 



where B = Vi2gW_ Having summed to all orders, the contributions of -^t+t in the N"LL contribution to V to 
obtain V//, we again find that V has a peculiar singularity structure. 



Appendix B: Massless Scalar Electrodynamics 



In this Appendix we show how the sums appearing in Eqs. ([551 [55)1 for Vll and Vnll can be evaluated in closed 
form by adapting the method of characteristics [TH [19] . This first entails defining 



'Wn+k {9{t),ot{t),t) = exp 



4 / 71 (ff(T),a(T)) dr 





pi+k (git), ait)) 



where 



dm 

dt 
da{t) 
dt 



(3l{g{t),m) (5(0) =ff) 
^^p^{g{t),a{t)) (a(0)=a) 



are characteristic functions. From Eqs. (fTUTlfTUgi) it follows that 
d 



dt 



"n+fc(5,a,i) = ^/3f(g,a)^ +^^(5r,a)^ +47i(g,a)^ w^+fc(g,a,i) 



Together, Eqs. jMl show that 

1 / 5 d 

Wn+i{g,a,t) = — f /3|(^,a)— + /?^(5,a)— +47i(g,a) ) u;""^(g,a,i) 

We now define 

oo 

VLLit) = Y,wl+, (g(i),a(i),i)ZV , 

where 

^ = loR, _2, . 



n=0 



(107) 

(108) 
(109) 

(110) 

(111) 

(112) 

(113) 



^Solutions to Eqs. I I108I [TOOt appear in Ref. [l]. They are easily obtained as b"j^ = 0. 
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with 



dp.{t) 
dt 



fi{t) , ^(0) = 



From Eqs. ([551 [TUTHTUgi fTTWlH) it follows that 

VLLit = 0) = Vll ■ 
We see that Eqs. pTUl [TTT|) lead to 



so that Eq. (|112p becomes 



n=0 



a U + 



L 



Furthermore, Eqs. pUTHTUgi [TTC)) reduce Eq. pT7)l to 



This coincides with the result appearing in Ref. [5]. 

Having found this closed form expression for the leading-log contribution to V{(f>), we turn to the next- 
log contribution of Eq. ([55]) . The first step is to define 

oo 

V^LLit) - <+2 (9{t), a{t), t) ^'V' • 

n=0 



(114) 
(115) 
(116) 

(117) 

(118) 
■to-leading 

(119) 



We now note that Eqs. [TTUl [TTT1) together show that 



^^71-1-2 



2n 



dt 



d 



d 



'■"^g^^^ da 



(120) 



where D{t) corresponds to the differential operator acting upon wJJ ^ . 
Iterating Eq. (|120p . we obtain 



1 

2n 



1 



dt \2{n- 1) \dt 
Repeating this n times and using Eq. (|116p we find 
1 



"+2 " 2"7i! 
The identity 



d" n d"-^ , , d''-"^ , ,d n -^""^ , . d^ „ 



(it" ^ 



dt 



dt"-3 



dt2 



D{t) 



dt- 



d" 
dt 



converts Eq. (|122p to the form 



2"n! 



|;»5 + ^(6»!)-5 



d" 



d"w? 



df^ 



where 



dt 



D^D 



(121) 



(122) 



(123) 



(124) 



(125) 
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More explicitly, we have 



(31 [5(r), a(r)] + . . . + 2 (71 [g(r), a(r)])^ 



i?(t)u;°(.9(t),a(t),i) = /3f[g(0,«W] 




2(7i[5W,«(t)]r 



(.g(t),a(t),t) 



(126) 



/3f [5(r),a(r)] 



2(71 [9{T),o^{r)]Y 



dT 



dt 



w1 {g{t),a{t),t) 



(127) 



Eqs. (|126p and (|127p ensure consistency between Eqs. (|122p and (|124p . In Eqs. (|126p and (|127p . 5 and a are evaluated 
at t when appearing in the arguments of . Derivatives with respect to g and a have these functions evaluated at 
the scale t. In Eq. (|127p , derivatives with respect to t acting on 'Wi{g{t), a{t), t) do so prior to functional derivatives 
d/dg{t), d/da{t)-, the last step in Eq. p27p is the integral over t. 
Eqs. (Hini [nil) lead to 



VNLLit) = (^g{t + a{t + ^),t+^^+D(^t + ^^ u;? (^g{t + a{t + ^),t+ ^) 



D{t)w°, ( g{t + ^),a{t + ^),t + ^ 



(128) 



It is evident that 

VNLLit = 0) VnlL 

and so Vnll of Eq. 
Vnll = 



^^(4),4),f)+5(f).o(,(f ),.(§),§ 



(129) 



(130) 



(Since D'it) = D{t), we can set D(0) = 0.) 

As a result, we see that if /3|, /Jf , /3f , 71, 72,_Pi and are known, then Vnll is fully determined. It is seen 
that this approach can also be used to find closed form expressions for V/v^ll(p > 2). 
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